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The generalized maximum entropy method (GMEM) is a special modification of the standard maximum 
entropy method (MEM) which seeks solutions in the space of complex functions. In this work a reduced 
version of the GMEM intended for reconstructing real images with positive and negative values is used. 
As compared with the standard MEM, intended for the reconstruction of only non-negative images, the 
GMEM allows us to obtain higher-quality images with a much lower level of nonlinear distortions caused 
by errors in the data. Here, we present the results of the GMEM imaging of 36 selected extragalactic 
radio sources with a resolution of 0.3-0.5 mas on astrometric and geodetic VLBI observations at 8.2 
GHz, obtained with a global array in the period from 1994-1996. In VLBI mapping practice this is the 
first experience of imaging with such a high resolution using maximum entropy technique. A differential 
maximum entropy method intended for increasing the dynamic range of images is demonstrated on the 
radio source 0059+581. In the case of unreliable 'closure' phases, completely 'phaseless' methods of 
mapping are recommended. Maps of two sources 0615+820 and 0642+214 are obtained using one such 
method. 

KEY WORDS VLBI imaging, generalized maximum entropy method, compact extragalactic radio 
sources 



1 INTRODUCTION 



There are two main deconvolution algorithms in Very Long Baseline Interferometry (VLBI) 
imaging: CLEAN and maximum entropy method (MEM). The maximum entropy method is 
well known among physicists as a very powerful tool for the nonlinear regularization of incorrect 
tasks. But in spite of its high super resolution effect, MEM is less popular in VLBI. There 
are two main reasons: the high computational complexity and the nonlinear image distortions 
caused by errors in data (visibility function). But the situation can be change fundamentally 
by introducing into VLBI imaging an effective modification of the MEM called the generalized 
maximum entropy method (GMEM) which generally operates in the space of complex functions 
(Bajkova, 1990, 1993). Searching for solutions of complex functions in space allows us firstly 
to factor multidimensional MEM algorithms into series of simpler one-dimensional algorithms, 
which considerably decreases the computational complexity (Frieden and Bajkova, 1994) and, 
secondly, considerably decreases the nonlinear distortions of images (Bajkova, 1995). 
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Thus, due to its generalized form MEM can become preferable in VLBI to CLEAN, especially 
for imaging sources with a complicated extended structure (we know, that MEM produces 
maximally smooth images, but CLEAN sharpens them). 

In this paper we present the first experience of applying the GMEM to mapping compact 
extragalactic sources with high resolution (0.3-0.5 mas). We used VLBI astrometric and geodetic 
observations (NEOS-A program) obtained with a global array at 8.2 GHz. Data were received 
from the Goddard Space Flight Center. Source imaging from astrometric and geodetic VLBI 
observations is of great interest both for astrometry, geodesy (geodynamics taking into account 
extended source structure on mas scales to more accurate determination of celestial/terrestrial 
coordinates) and astrophysics (investigation of the short-period structure variability of compact 
extragalactic radio sources). 

For mapping we used the well-known 'CalTech VLBI Program' package into which we intro- 
duced the GMEM procedure as a deconvolution operation in the selfcalibration loop. 

2 MAPPING COMPACT EXTRAGALACTIC RADIO SOURCES 



For VLBI source mapping a reduced version of the GMEM intended for the reconstruction 
of real signals with positive and negative values was used. The selection of the GMEM instead 
of the MEM as a deconvolution operation was principally dictated by low signal-to-noise ratio 
typical of astrometric and geodetic VLBI observations. 

A detailed description of the practical GMEM algorithm is given in Bajkova (1998a). The 
size of maps was chosen as 256x256 pixels, with an interval between two pixels of 0.1 mas and 
an image domain size between 60x60 and 100x100 pixels, depending in each case on structure 
extent and data quality. Parameter 'a' (Bajkova, 1990, 1998a) in the GMEM algorithm, which 
is responsible for separating 'positive' and 'negative' solutions, was chosen equal to 10 8 . A 
numerical solution of the entropic functional was realized using steepest-descent method. To 
ensure the high stability of the numerical algorithm to noise, an additional regularizing term 
was added to the principal entropic functional (Bajkova, 1998a). 

Images of 36 selected compact extragalactic radio sources are shown in Figure 1. Character- 
istics of the sources are given in Table 1. Parameters of the maps (date of observations, list of 
VLBI stations, flux densities and agreement factors) are given in Table 2. The presented images 
were obtained by convolution of entropic solutions by a circular gaussian beam with a circle 
diameter on 0.2 mas. This size is enough to smooth over image samples distant which are 0.1 
mas from each other (the beam circles are shown in bottom left corner of the maps). 

The GMEM images can be compared with the corresponding CLEAN images published 
in Bajkova et al. (1996a). The comparison shows, that the GMEM maps are obtained with 
better agreement factors (nearly in 2.2 times) and are less 'lumpy'. In this connection the total 
flux increased on the average by 1.036 times, i.e. practically did not change. But the peak 
flux decreased by 1.83 times. Taking into account this and the fact that MEM and GMEM 
principally seek a solution in space of smooth gaussian functions, it is possible to conclude that 
the new GMEM maps are smoother than the CLEAN maps. 
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Improvement of agreement factors in case of GMEM shows that the space of smooth functions 
for searching for source distributions is more adequate than the space of S— functions in which the 
CLEAN solutions are sought, i.e. the considered sources reveal rather extended structure at the 
given resolution. Note that CLEAN always gives disconnected solution as a combination of 5— 
functions (in addition, CLEAN sharpens bright components and damps weak ones), although 
CLEAN maps usually seem smooth. This is because a CLEAN map is a convolution of the 
CLEAN solution (model) by a smooth gaussian 'clean' beam. And, obviously, in general a 
CLEAN map is not obliged to agree with the data to such an extent as the corresponding 
CLEAN model. But when we use the GMEM we at once obtain a smooth solution which is in 
good agreement with the data. 
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Table 2. Parameters of maps 



Source Date Stations Flux [Jy] Agr. factor 
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Abbreviations: A - Algopark; BR - VLBA-BR; C - Crimea; D65 - DSS65; F - Fortleza; G - 
Gilcreek; H - Hartrao; H26 - Hobart26; K - Kokee; MA - Matera; ME - Medicina; MK 
VLBA-MK; N - Noto; N20 - NRAO20; N85 - NRA085; NL - VLBA-NL; NY - NyAlesund20; O 
-Onsala; SC - VLBA-SC; W - Wettzell. 



If we obtain GMEM solution with a worse agreement factor than the CLEAN one, we 
should conclude that in this case GMEM is less adequate than CLEAN and the source can be 
considered to be point-like source consisting of some 5— functions. In our case, such a source is 
0642+214 (in fact, its CLEAN solution has better agreement factor than the GMEM solution 
(by approximately two times)). 
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On the basis of careful analysis of both the CLEAN and GMEM maps we can conclude that 
the GMEM gave more accurate estimates of source distributions than the CLEAN map. In some 
of the sources (0400+258, 0458-020, 0917+624, 1606+106, 1823+568, 2230+114) the extended 
structure details are revealed to a fuller extent. All the images are free from the artifacts inherent 
in CLEAN. 

To increase the dynamic range of the maps we used a differential GMEM (see Section 3) . A 
more detailed comparison of CLEAN and GMEM images is given in Bajkova (1999). 

It is also necessary to note that the maps of two sources with unreliable 'closure' phases 
0615+820 and 0642+214 were obtained using completely 'phaseless' mapping (see section 4). 

3 A DIFFERENTIAL METHOD OF MAXIMUM ENTROPY 



At present, the idea of differential mapping is realized most comprehensively in the 'DIFMAP' 
program package at the CalTech, where CLEAN method is used as the deconvolution operation. 

The method of differential mapping is based on the fundamental linear property of linearity 
of the Fourier transformation. According to this method, the bright components of the source, 
which were reconstructed at the first stage, are subtracted from the initial visibility function, 
the subsequent reconstruction is performed using the residual visibility function, and the recon- 
struction results are added at the final stage. 

In the case of CLEAN, such a mapping method mainly influences the convergence rate. 
However, if the MEM, which has pronounced nonlinear properties, is used as the deconvolution 
operation, the principle of differential mapping can improve the reconstruction quality, particu- 
larly when the source has bright compact components against the background of a sufficiently 
weak extended base. 

What is the cause of the reconstruction quality improvement? The fact is that after the 
subtraction of bright components, which were reconstructed at the first stage using the MEM, 
from the initial visibility function we obtain the residual visibility function, in which the pro- 
portion of the weak extended component becomes larger. Therefore, we artificially decrease the 
dynamic interval of the map, which corresponds to the residual visibility function, thus simplify 
the image reconstruction at the second stage. 

Above we discussed the advantages of the GMEM for processing real non-negative images. 
With respect to differential mapping, another advantage of GMEM originated from the fact that 
an image with negative values can correspond to the residual visibility function after subtraction 
of the bright component that is reconstructed at the first stage of the algorithm from the initial 
visibility function. This is possible if the bright component was reconstructed at the first stage 
with an overestimated amplitude, which is quite typical of any nonlinear method. To rule out 
unwanted image distortions when a non-negative solution is sought using the data assuming 
the presence of negative components, we should use the generalized method rather than the 
standard one. 

Simulation results of the differential GMEM are shown in Bajkova (1998b). Here we demon- 
strate potential of differential GMEM for mapping the radio source 0059+581, which is actively 
used as a reference source in the astrometric and geodetic VLBI programs. This source is of 
interest because it shows fast variation of both the total flux and the structure, which consists of 
a compact core and rather weak extended elements (jet). The maps that were constructed using 
GMEM for a number of dates between June 1994 and December 1995 are shown in Bajkova 
(1998b). Between the middle and the end of 1994 we observed an almost linear decrease in the 
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total flux from ~ 4 Jy to 1.5 Jy. As is obvious from Bajkova (1998b), we failed to obtain a 
map with a sufficient dynamic range detect an extended element from June 26, 1994, when the 
flux was maximal over the interval in question. Using the standard GMEM, we only detected 
extended elements from October 4, 1994, when their share with respect to the total flux became 
sufficient to be detected by the selected reconstruction method. In Figure 2, for comparison, we 
show maps, obtained using both the standard and differential GMEM for three dates on which 
the total flux was sufficiently large and the flux share corresponding to extended elements was 
insignificant. As is obvious from the figure, the use of differential mapping allowed us to increase 
the dynamic range of the maps such that the extended structures became very distinguished. 



4 'PHASELESS' MAPPING 



A principle of 'phaseless' mapping used for imaging two sources 0615+820 and 0642+214 
(see Figure 1) is based on reconstruction of the visibility function amplitude over the whole 
UV-plane. The image is then reconstructed using only the amplitude of the visibility function 
(the well-known 'phase problem'). For reconstruction of the visibility amplitude GMEM has 
also been used in the self-calibration loop setting all the 'closure' phases to zero. For solving the 
'phase problem' an effective modification of the well-known Fienup algorithm was used (Bajkova, 
1996b). 'Phaseless' mapping is recommended in the case of bad 'closure' phases. 

5 SUMMARY 



This work is devoted to the first experience of the high resolution imaging of compact extra- 
galactic sources using generalized maximum entropy method proposed by the author in earlier 
papers (Bajkova, 1990, 1993). Images of a sample of 36 selected compact extragalactic radio 
sources obtained with a resolution of 0.3-0.5 mas on astrometric and geodetic VLBI observations 
at 8.2 GHz with a global array (NEOS-A program) between 1994 and 1996 are presented. Com- 
parison with the corresponding CLEAN images shows that the images obtained are smoother 
and characterized by better agreement factors, revealing rather complicated structures to a fuller 
extent. 

The choice of reconstruction method as the deconvolution procedure in selfcalibration loop 
depends on the signal-to-noise ratio and a priori information about the source's structure. In the 
case of point-like structures, which are the majority of sources observed in geodetic programs, the 
traditional CLEAN is the preferable procedure because it was invented especially for such sources 
and is very fast. In the case of complicated structures with strongly extended components, the 
nonlinear information methods can be preferable. Moreover, in the case of a low signal-to- 
noise ratio, which is typical for astrometric and geodetic VLBI observations, the generalized 
information methods for seeking a solution in the space of functions with positive and negative 
values are preferable to the standard ones which deal only with real non-negative distributions, 
because they ensure a much lower level of nonlinear distortions. 

To increase the dynamic range of images we used a differential generalized maximum entropy 
method. Compared with the traditional methods, the nonlinear methods of differential mapping 
using the maximum entropy method as the deconvolution operation allow us to substantially 
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improve the reconstruction quality of source images containing bright point components against 
a background of a weak base. The maps obtained using the differential method of maximum 
entropy are characterized by a higher dynamic range. To eliminate possible nonlinear distortions 
in the case of differential mapping, the generalized method of maximum entropy is preferred to 
the standard method. 

In the case of unreliable 'closure' phases we recommend using completely 'phaseless' methods 
of mapping. 
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Figure 1 Images of extragalactic sources obtained by GMEM. Contour levels are 1, 2, 4, 10, 20, 30, 40, 
50, 60, 70, 80, 90, 99 % of peak intensity. 
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Figure 1 (Continued). 
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Figure 1 (Continued). 
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Figure 2 Images of source 0059+581 obtained by the standard GMEM (top) and the differential GMEM 
(bottom) for three dates. Contour levels are 1, 2, 4, 10, 20, 30, 40, 50, 60, 70, 80, 90, 99 % of peak 
intensity. 



